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Lattice Gauge Theory enables an ab initio study of the low-energy properties of 
Quantum Chromodynamics, the theory of the strong interaction. I begin these 
lectures by presenting the lattice formulation of QCD, and then outline the bench- 
mark calculation of lattice QCD, the light-hadron spectrum. I then proceed to 
explore the predictive power of lattice QCD, in particular as it pertains to hadronic 
physics. I will discuss the spectrum of glueballs, exotics and excited states, before 
investigating the study of form factors and structure functions. I will conclude 
by showing how lattice QCD can be used to study multi-hadron systems, and in 
particular provide insight into the nucleon-nucleon interaction. 



1 Lattice QCD: the Basics 

1.1 Introduction 

The fundamental forces of nature can by characterised by the strength of the 
interaction: gravity, the weak interaction, responsible for /3-decay, the electro- 
magnetic interaction, and finally the strong nuclear force. All but the weakest 
of these, gravity, are incorporated in the Standard Model of particle interac- 
tions. 

The Standard Model describes interactions through gauge theories, char- 
acterised by a local symmetry, or gauge invariance. The simplest is the electro- 
magnetic interaction, with the Abelian symmetry of the gauge group U(l). The 
model of Glashow, Weinberg and Salam unified the electromagnetic and weak 
interactions through the a broken symmetry group SU(2)(g>U(l). The strong in- 
teraction is associated with the unbroken non- Abelian symmetry group SU(3), 
and is accorded the name Quantum Chromodynamics (QCD). 

The strength of the electromagnetic interaction is characterised by the 
dimensionless fine-structure constant a e ~ 1/137. A very powerful calcula- 
tional technique is to expand as a series in a e - perturbation theory. QCD 
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is characterised by a strong coupling constant a s ~ C(l). QCD, however, 
is asymptotically free, with an effective running coupling a s (Q 2 ) decreasing 
logarithmically with increasing Q 2 . Thus processes with an energy scale large 
compared with the natural scale of the strong interaction, of the order of the 
proton mass, are often amenable to the techniques of perturbation theory. 

At energy scales of the order of the proton mass, perturbation theory fails. 
Yet a quantitative understanding of QCD is crucial both for the study of the 
strong interaction, and for the study of the other forces which are masked 
by the strong interaction. In this energy regime, we can either employ low- 
energy effective models of QCD, or seek some way of performing a quantitative 
calculation directly within QCD. Lattice QCD is the only means we have of 
performing such an ab initio calculation. 

Before proceeding to a description of lattice QCD, it is useful to make a 
comparison between the properties of QCD and of QED: 



The gluon self-coupling reflects the non-Abelian, and highly non-linear, nature 
of QCD. Where are the quarks? They are bound into the colour-singlet hadrons 
we observe in nature. Lattice gauge theory provides the means to relate the 
quark degrees of freedom with the observed hadronic degrees of freedom. 

1.2 Lattice Gauge Theory 

Lattice gauge theory was proposed by Ken Wilson in 19740 Because of the 
gluon self-coupling, we have a sensible pure-gauge theory of interacting gluons, 
even without quark, or matter, fields. We will consider this theory first. 

We begin by formulating QCD in Euclidean space, which we accomplish 
by a Wick rotation from Minkowski space, 
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Coupling to 
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(1) 



The gauge fields are defined through 




(2) 



[T a ,T b ] = iffT c 



(3) 
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Figure 1: A schematic of a lattice showing the association of the SU(3) matrices U^ix) with 
the links of the lattice. 



Tr T a T b = -S ab . (4) 

We now introduce the field-strength tensor 

F a ^ = d^Al-d v Al+gf abc AlAl, (5) 
in terms of which the Euclidean continuum action is 

S=\J 'd 4 xsF^ v . (6) 

As we will see later, the crucial property of Euclidean space QCD for the 
formulation of lattice gauge theories is that the action is real. Gauge invariance 
is manifest through invariance under the transformation 

A„(x) -> A(x)A^(x)A- 1 (x) - -(cVK^A- 1 ^). (7) 

We proceed to the lattice formulation of QCD by replacing a finite region 
of continuum space-time by a discrete four-dimensional lattice, or grid, of 
points. The gluon degrees of freedom are represented by SU(3) matrices U^{x) 
associated with the links connecting the grid points, as shown in Figure m 
We work with the elements of the group, rather than elements of the algebra, 
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and the SU(3) matrices U^(x) are related to the usual continuum gauge fields 
through 



Ufj,{x) = expiga / dt A fl (x + tafi), 



(8) 



where g is the coupling constant, and a the lattice spacing. Under a gauge 
transformation A(x), the link variables transform as 



U^x) -» A(x)C/ M (x + A)A- 1 (x), 



(9) 



in analogy with Eq. (?]. Wilson'sJorm of the lattice gauge action is constructed 
from the elementary plaquettesu 

Un^{x) = Up{x)U v {x + jjL)Ul{x + v)U$(x). (10) 

The plaqucttes are clearly gauge invariant, and the action is then written 



X (A>V 



l-l-STrUn^x) 



= -#EE KTr ^ en) 



N, 



X Ll>V 

where we have ignored the constant term, and introduced 
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with, for QCD, iV c = 3. It is straightforward to show that the Wilson lattice 
gauge action is related to the continuum counterpart, Eq. (^), by 

S G = ^Jd 4 xF^F^+0(a 2 ), (12) 

so that the lattice gauge action has C(a 2 ) discretisation errors. 

1.3 Observables and Lattice Gauge Simulations 

Within lattice gauge theory, the expectation value of an observable O is given 
by the path integral 



where 



(°) = |/ VUO(U)e 
VU = Y[dU^(x) 



Sg(U) 



(13) 
(14) 
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and Z is the generating functional 

Z = f VUe- Sa( - u) ;VU = Y[dUf,{x). (15) 

Before proceeding further, we need to define what we mean by the integration 
over a group variable dU . We do this through the Haar measure, which for a 
compact group is the unique measure having the following properties: 
1. 

/ dUf(U)= [ dUf(VU)= [ dUf(UV) WeG. 
Jg Jg Jg 

2. 

/ dU = l. 
Jg 

This choice of measure respects gauge invariance. Note that, because we are 
employing the compact variables, U ll {x) 1 rather than the elements of the alge- 
bra, we do not need to fix the gauge, and indeed in most circumstances we do 
not do so. However, there are cases where working in a fixed gauge is useful, 
most notably in lattice perturbation theory, where gauge fixing is essential, in 
the definition of hadronic wave functions, where it is often useful to work in 
Coulomb gauge, and most directly in the study of the fundamental gluon and 
quark Green functions of the theory. 

On a finite lattice, the calculation of observables is equivalent to the evalu- 
ation of a very high, though finite, dimensional integral. In principle, we could 
estimate this integral by evaluating the integrand at uniformly distributed 
points. This, however, would be hopelessly inefficient; the exponential be- 
haviour ensures that the integral is dominated by regions where the action is 
small. Instead we use importance sampling, and generate gauge fields with a 
probability distribution 

e~ SG(u \ (16) 

The interpretation of this exponential in terms of a probability distribution 
requires that the action be real, and hence the need to work in Euclidean 
space. The formulation follows that of many systems in statistical physics. 

1.4 Statistical and Systematic Uncertainties 
Statistical Uncertainties 

Observables in lattice QCD calculations arise from a Monte Carlo procedure, 
and thus have statistical uncertainties. Once we have reached thermalisation, 
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these uncertainties decrease as the square root of the number of configura- 
tions, providing successive configurations are sufficiently widely separated to 
be statistically independent. 

Systematic Uncertainties 

Of even greater delicacy than the statistical uncertainties are the system- 
atic uncertainties that enter our computations. These arise from a variety 
of sources, including: 

• Finite Volume: Our box must be sufficiently large that finite volume 
effects are under control. For light hadron spectroscopy, box sizes of at 
least 2 fm are necessary to ensure that the hadron is not "squeezed" , but 
for excited states even larger volumes may be required. In addition, the 
requirement that the spatial extent of the lattice be large compared with 
the correlation length, set by the pseudoscalar mass, sets a still more 
stringent constraint at the physical pion mass. 

• Discretisation Effects: Increasing the inverse coupling (3 corresponds to 
progressing to weaker coupling, and hence smaller lattice spacing a. We 
must ensure that (3 is sufficiently large that the scale-breaking discreti- 
sation errors arc under control, and in practice we perform calculations 
at several values of a and extrapolate to the limit a = 0. 

We will encounter several other potential sources of systematic errors when we 
discuss the inclusion of the quarks. 

1.5 Including the Quarks 

The full generating functional for lattice QCD with a single flavour of quark is 

Z = J TO^^e- 5 ^^^-^^^*^^^, (17) 

where M(x,y, U) is the fermion matrix which, in its "naive" form, is 

M(x,y,U) =m5 Xi y + ^^7^ (£^( X )<W+A ~ u t( x - £)<W-a) ( 18 ) 

with m the quark mass. Because the fermion fields are represented by Grass- 
man variables, we can integrate out the fermion degrees of freedom, to obtain 

Z = [ VU dctM(U)e~ Sc{u) . (19) 
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The determinant fluctuates rapidly between configurations. Thus it is not suf- 
ficient for a Monte-Carlo procedure to generate configurations with probability 
e~ s °( u \ and only include the determinant in the calculation of observables; 
det M(U) must be included in the measure. Furthermore, whilst multiplica- 
tion by the fermion matrix M involves only nearest neighbour communication, 
the evaluation of det M(U) is essentially a global operation. Thus detM(J7) 
must be re-evaluated every time we update even a single value link variable 
U^(x). The most efficient algorithms for the simulation of QCD with dynam- 
ical quarks, such as the Hybrid Monte Carlo algorithmp involve a non-local 
updating procedure. Nevertheless, calculations with dynamical quarks are at 
least 1000 times as expensive as pure gauge calculations. 

The computational overhead of completely including the quark degrees of 
freedom has encouraged many calculations to be performed in the quenched 
approximation to QCD, in which we set 



in Eq. ([19]). This corresponds to suppressing the contribution of closed quark 
loops in the path integral. There are two justifications for this seemly radical 
approximation. Firstly, the phenomenological observation that the neglect 
of the quark loops corresponds to the neglect of OZFsuppressed processes. 
Secondly, the quenched approximation emerges in the large N c limit of QCD, 
where N c is the number of colours. 



1.6 Fermion doubling and chiral symmetry 

Unfortunately, the lattice formulation of fermions presents a further challenge. 
To illustrate how this arises, let us consider the momentum-space fermion 
propagator 



In the massless limit, the propagator has a pole not only at p^ = 0, but also 
at the edges of the Brillouin zone p^ = n/a. Thus in four dimensions, we 
have a theory with 2 4 non-interacting, equal mass fermions. This situation is 
a consequence of the Dirac equation being first order. Historically, there have 
been two solutions to this problem. 



detM = 1, 
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Wilson Fermions 



Wilson proposed the addition of a second derivative term, or momentum- 
dependent mass term, to the action: 

S^ = ^{(m + 4r)#r)#r)- 

X 

X - ^2 [${x){r - lll )U ll {x)^{x + A) + ${x + A)(r + iJUfaMx)] j .(21) 
In the continuum limit, we find 

/orD 2 
d 4 x$(x)(D + m —)^{ X ) + C{a 2 ). (22) 

The addition of the second derivatives lifts the mass of the unwanted doublcrs, 
but at a price. We have added 0(a) discretisation errors to the fermion matrix, 
and furthermore the additional term explicitly breaks chiral symmetry at any 
non-zero value of the lattice spacing, though it is important to remember that 
chiral symmetry is restored in the continuum limit. The breaking of chiral 
symmetry has the unfortunate consequence that the fermion masses are sub- 
ject to an additive mass renormalisation. In simulations, it is conventional to 
reparameterise the fermion matrix as 

M x , y = (m + 4r) {5 x . y l - n x "hopping term"} (23) 

where the hopping parameter 

k = — — ^ r (24) 

2(4r + to) v ; 

is now a tunable parameter, reflecting the additive quark-mass renormalisation. 
The critical value of the hopping parameter, K cr ;t, is that value for which the 
pion mass vanishes. 



Kogut-Susskind Fermions 

The second approach, due to Kogut and Susskind, regards the complete loss 
of chiral symmetry at finite a as too great a sacrifice, and therefore aims to 
preserve a remnant of chiral symmetry whilst reducing the flavour-doubling 
problem. This formulation leads to four flavours of quark, but the different 
flavour and spin components are assembled from fields at the sixteen corners of 
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a 2 4 hypercube. While the Kogut-Susskind formulation has been extensively 
employed in the calculation of quantities where chiral symmetry is crucial, the 
problematic flavour identification means that I will concentrate on the results 
using the Wilson formulation in these lectures. 

Chiral Fermions and the Ginsparg- Wilson Relation 

Is it possible to construct a formulation that does indeed possess a symmetry 
analogous to chiral symmetry at a finite lattice spacing, whilst admitting the 
correct spectrum of quark states? Let us list four properties desired of the free 
fermion action, which we write in the form 

S F = a A ^2'ii>{x)D{x 1 y)'4){y) 

X 

1. D(x, y) is local 

2. Far below the cut-off, D{p) ~ ij^Pp, + 0{p 2 ). 

3. D(p) is invertible at all non-zero momenta. 

4. l5 D + D Tj = 0. 

The last two requirements require explanation; ^ demands that the only poles 
occur at zero momentum, and hence that there be no doublers, whilst ^ is 
just a statement of chiral symmetry. The famous Nielsen-Ninomiya "no-go" 
theorem Q states that it is not possible to find a Dirac operator that allows 
all four requirements to be satisfied simultaneously, and since the first two 
requirements were deemed sacrosanct either chiral symmetry had to be broken, 
or flavour-doubling had to be accepted. 

The recent revolution in our understanding of_chiral fermions are related 
to the rediscovery of the Ginsparg- Wilson relation □ 

lb D + D l5 = aD l5 D, (25) 

as a replacement for requirement ^. At non-zero distances, the Dirac operator 
does indeed commute with 75 , and a symmetry, reducing to chiral symmetry in 
the continuum limit, is preserved even at non-zero values of the lattice spacing. 
The problem is finding a formulation of the Dirac operator that does indeed 
satisfy the relation of Eq. |25|. 

Recently, two formulations satisfying this relation have been discovered. 
In the case of Domain Wall fermionsM (DWF), an auxiliary fifth dimension, 
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Figure 2: The residual pjen mass as a function of the extent L 3 of the fifth dimension 
in the DWF formulationQ The calculation is performed in the quenched approximation to 
QCD, and Wilson and Iwasaki refer to the standard Wilson gauge action, and an improved 
gauge action respectively. 

with coordinate s and extent N s , is introduced. The action is essentially a 
five-dimensional Wilson fermion action 

5° WF = Yl Hx)(D(x )y )S s , s ,+D 5 (s, S ')5 x , v )^(y), (26) 

where D(x,y) is the usual Wilson-Dirac operator introduced in Eq. ||l], but 
with a negative mass term M. The operator in the fifth dimension, _D 5 (s,s'), 
couples the boundaries through a parameter —to, which is proportional to the 
usual four-dimensional quark mass; note that no gauge links are introduced in 
the fifth dimension. The chiral limit corresponds to L s — > oo, and then to — > 0; 
following the former of these limits, the quark mass is only multiplicatively 
renormalised. 

A crucial issue is how small a value of L s is sufficient to maintain good 
chiral properties whilst minimising the computational cost. For the case of 
hadronic physics, perhaps more important than having good chiral properties 
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is being able to perform simulations at sufficiently small quark mass for the 
pion cloud to emerge. Fig. ^ shows the residual pion mass in the chiral limit 
as a function of function of L s u 

The second, approach is the Overlap formalism, introduced by Narayanan 
and NeubergerEra Here the overlap-Dirac operator satisfying the property [| is 

^ = ^ + M+(l-M)75 f ( "> (27) 

where H(m) is the Hcrmitian Wilson-Dirac fermion operator, with negative 
mass m, defined by H = 75-D where D is the usual Wilson-Dirac operator. 
The parameter fi is related to the physical quark mass. In this case, the extra 
computational cost comes not from computing in five dimensions, but rather 
from evaluating the step function 

e(H) = -£=. (28) 

The relative computational .overheads of the two implementations is the sub- 
ject of intense investigational but in any case the overhead compared to the 
standard Wilson fermion action is considerable. Whilst this overhead is justi- 
fied for chiral gauge theories, the situation in the case of hadronic physics is 
less clear; perhaps there are more efficient ways of approaching physical values 
of the light-quark masses. 



1.7 Improvement 

The addition of the Wilson term to the fermion action has introduced 0(a) 
discretisation errors; in contrast the gauge action has only 0(a 2 ) discretisation 
errors. Thus there has been an emphasis on reducing the errors in the fermion 
sector through the addition of higher Jimcnsional terms to the action, the 
improvement programme of SymanzikEU In the case of the Wilson fermion 
action, the leading 0(a) errors can be removed through the addition of a 
single dimension-five operator, the magnetic moment, or clover term, proposed 
by Sheikholcslami and WohlertEj 

Sf = S f - J2 ^)F^{x)a^(x). (29) 

The name "clover" is clear from the natural lattice discretisation of illus- 
trated in Figure ^. 
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Figure 3: The lattice discretisation of the field-strength F^lx) in terms of the plaquettes 
with corners at x 



Using tree-level perturbation theory, the clover coefficient csw is unity, and 
the discretisation errors on hadron masses and, with an appropriate discreti- 
sation of operators, on-shell matrix elements are formally (D(ag 2 ).£3 UKQCD 
performed an extensive investigation of the hadron spectrum and hadronic ma- 
trix elements using this value of Csw ; and the discretisatiqn.errors, particularly 
for systems containing heavy quarks, can be substantial. Efl 

More recently, two prescriptions for determining csw have been proposed 
with the aim of reducing discretisation errors still further. In the first, the 
clover coefficient is constrained to its mean-field-improved, or tadpole, valucta 

c sw = TAD = i (30) 
u Q 

where 

uo = <|TH7 n > (31) 

is an estimate of the mean- value of the link variable [/„. Though formally 
the discretisation errors remain 0(ag 2 ), this prescription recognises the poor 
behaviour of naive lattice perturbation theory arising from the "tadpole" con- 
tributions, and attempts to resum the dominant higher-order contributions 
through the use of a more physical expansion parameter. 

The second prescription t3l!3 determines csw non-perturbatively in such a 
way as to remove all calO(a) discretisation errors from hadron masses, and, 
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with an appropriate choice of operators, from all on-shcll matrix elements El 

NP l- 0.656. 9o 2 -0.152. 9 4-Q.054ff 6 2 
c sw = NP = 2 , g < 1, (32) 



where gg = 6//?. 

We end this section by remarking that the chiral-fermion formulations are 
already C(a)-improved; the 0(a) discretisation are introduced through the 
chiral-symmetry-breaking Wilson term. 



2 The light-hadron spectrum 

The calculation of the spectrum of hadrons containing the light quarks (u,d,s) 
is the benchmark calculation of lattice QCD; we know many of the results! It 
also provides a useful theatre for discussion of some of the issues I raised in 
the introduction. However, let us begin with one of the simplest observations 
we can make in lattice QCD, that of the linear confining potential. 



2.1 The Static Quark Potential and Quark Confinement 

The simplest observable we can obtain from a lattice simulation is the potential 
between two (infinitely heavy) static quarks. We construct the Wilson loops 

W(R, T) = (Tr U(C(R, T))) (33) 

where U(C(R,T)) is the product of gauge links around a R x T space-time 
loop. At large times, we can extract the potential V(R) between two static 
quarks Q at separation R, using 

W(R,T) T ~°° e~ TV{ ~ R \ (34) 

An area-law decay of the large Wilson loops is characteristic of a linear, confin- 
ing potential in QCD, giving rise to a constant force with increasing separations 
r. We parameterise this force F(r) by 

F( r y = lM + ^(^-lj (35) 

where ro ~ 0.5 fm is a phenomenological parameter which we use to deter- 
mine the scale. We show this force for the pure-gauge theory, corresponding 
to the quenched approximation, in Figure [s]. At large separations we see the 
constant force indicative of confinement. But there is a further important ob- 
servation. The results from all three calculations lie on a single curve, even 



13 



C(R,T) 





Figure 4: The construction of the R X T wilson loop in a space-time plane. 



though the calculations span a factor of two in lattice spacing a, from about 
0.05 fm to 0.1 fm. 

Do we expect the same picture of a rising, linear static-quark potential in 
full QCD, with dynamical quarks? As the two heavy quarks are separated, 
the energy stored in the string increases. Eventually, the string can "break" to 
form a quark-antiquark pair, a process that does not occur in the pure-gauge 
theory. This should lead to a flattening of the potential at some distance 
corresponding to an energy in the flux tube of 2Mb, twice the binding energy 
of static-heavy-light meson. In practice, there has been no clear observation 
of such a feature at zero temperature from the simple Wilson loop operator; 
the behaviour of the potential in full QCD from a calculation using an 0(a)- 
improved fermion action Eil is shown in Figure ^. 

As the string breaks, the QQ system crosses to a system of two heavy- 
light Qq mesons. The lack of observation of the string breaking from the 
Wilson loop is ascribed to the poor overlap of the Wilson loop operator with 
two such heavy-light mesons. The mixing between these two states has been 
successfully investigatedJn simpler systems E3oO, and recently a study has 
been made within QCDt3 Both the ground state and first excited state energies 
are extracted by a variational calculation using the Wilson loop and QqQq 
operators, as shown in Figure ^. We will study a related problem when we 
discuss the nucleon-nucleon interaction in Section m 
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Figure 6: The scaled and normalised potential as a function of r/ro, where ro is the Sommer 
scale disortssed earlier, as obtained on a 12 3 X 24 lattice, using O (a)-improved dynamical 
fcrmions.EH The expected region of string breaking is shown by the horizontal lines. 
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Figure 7: The ground and excited state-ienergies obtained from a variational calculation 
including the QQ and QqQa, operatorsES The horizontal line is 2Mg, the static binding 
energy shown also in Figure pi 
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2.2 Spectrum recipe 

The elemental building blocks of a spectrum calculation are the quark propa- 
gators 



The quark propagator to every point x on the lattice from a fixed source point 
y, or linear combination of source points, is obtained by inverting the fermion 
matrix for a fixed source vector. There are a variety of linear-solver methods 
used to accomplish this. 

In principle, the recipe for determining the mass of a ground-state hadron 
P is straightforward: 

1. Choose an interpolating operator O that has a good overlap with P 



and ideally a small overlap with other states having the same quantum 
numbers. 

2. Form the time-sliced correlation function 



This is expressed using a Wick expansion in terms of the quark propa- 



3. Insert a complete set of states between O and . The time-sliced sum 
puts the intermediate states at rest, and we find 



G*U X ,y) = (Omx)i?Jy)\0). 



(36) 



(0 | O | P) ± 0, 



gators of Eq. 



X P 



d 3 k 



jr(0\O(x,t)\P(k))(P(k)\OH0,0)\0) 



(2Tr) 3 2E(k) 




p 



4. 



Continue to Euclidean space t — » it, and we find 




(37) 
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At large times, the lightest state dominates the spectral sum in Eq. (J37]), and 
we can extract the ground state mass. 

However there are many considerations that complicate this picture. Firstly, 
the temporal separation between the hadrons must be sufficiently large that 
the ground state can be identified in Eq. (|37|); we aim to accomplish this by 
choosing an operator having a large overlap with the ground state relative to 
the excited states, and by fitting to several interpolating operators. Secondly, 
the correlation lengths in our calculation must be small compared with the 
size of the box in which we are working. This correlation length is simply 
the inverse of the pion mass, and we require m^L ~ 5, where L is the spatial 
extent of the lattice. Most, though not all, simulations to date have restricted 
consideration to quark masses in the region of the strange-quark mass. 

The benchmark calculation of the quenched light-hadron spectrum us- 
ing the unimproved fermion action has been performed by the CP-PACS 
collaborationE3 They perform their calculation on a variety of lattice sizes, 
to control finite size effects, and at a variety of lattice spacings, to enable an 
extrapolation to the continuum limit. 

Recently, there has been a similar calculation using the noiL-perturbatively 
improved clover fermion action by the UKQCD CollaborationEZl They also per- 
form an extrapolation to the continuum limit, but in their case the discretisa- 
tion uncertainties are 0(a) rather than C(a 2 ), and I will describe some of the 
details of this calculation. 

UKQCD generated 16 3 x 32 lattices at j3 = 5.7, 16 3 x 48 lattices at /3 = 6.0 
and 24 3 x 48 lattices at f3 = 6.2, corresponding to approximately a factor 
two span of lattice spacings but at roughly the same spatial volumes. In 
addition, a calculation was performed on a larger 32 3 x 64 lattice at j3 = 6.0 
to enable a study of finite-volume effects. Quark propagators were computed 
with the clover coefficient having both its tadpole-improved value csw = TAD 
(J3 = 5.7,6.0,6.2), and its non-perturbatively determined value csw = NP 
((3 = 6.0,6.2). 

Since the quark propagators were computed at values of the quark mass 
in the region of the strange quark mass, it is necessary to extrapolate in the 
quark mass to obtain results at the u- and d-quark masses, and interpolate to 
obtain results at the s-quark mass. The bare, or unrenormalised quark mass, 
is related to the hopping parameter k through 

m q = f ( - - —) . (38) 

The bare quark mass must be rescaled in the C(a)-improved theory so that 
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spectral quantities approach the continuum limit with O(a 2 )0 

rh q = mi(l + b m am q ), (39) 

where the perturbative one- loop value of b m is usedil To extrapolate the results 
to the physical quark masses, the following ansatz is employed: 



m 



,2 
PS 



= B(m qt i + rh qa (40) 

my = Ay + Cy(m q ^ + m q ^) (41) 

TOOct = A 0c t + Coct(m q ,i + m q ,2 + m q ^) (42) 

m Doc = ^4 Dcc + C Dec (mg,i + m g ,2 + m,, 3 ), (43) 

where the subscripts PS, V, Oct and Dec refer to the pseudoscalar meson, 
vector meson, octet baryons (E- and A-like) and decuplet baryons (A-like) 
respectively. The critical hopping parameter corresponds to the value of K for 
which tops vanishes. The quality of the chiral extrapolations of the vector, E 
and A masses for the csw = NP data is shown in Figure ||. 

The csw = TAD data has in principle a remnant 0{a) discretisation error, 
whilst the csw = NP data is fully C(a 2 )-improved. In the continuum extrap- 
olations, UKQCD performs a simultaneous fit to both the NP and TAD data. 
In order to investigate the approach to the continuum limit, it is not neces- 
sary to study the chirally extrapolated values. Indeed, a clear demonstration 
of the efficacy of improvement can be seen by looking at the lattice-spacing 
dependence of hadron masses at a fixed m^jm p ratiopj shown in Figure^. 

The final UKQCD results for the quenched light-hadron spectrum, to- 
gether with the CP-PACS results using the standard Wilson fermion action, 
are shown in Figure The different UKQCD plotting symbols correspond 
to determining the lattice spacing by requiring that cither K* or N have its 
physical value. Both calculations support the assertion that the quenched 
light-hadron spectrum agrees with experiment at the 10% level. Now that we 
have established the reliability with which we can compute the known hadron 
masses, it is time to investigate the predictive power of lattice QCD. We begin 
by considering the glueball spectrum. 



2.3 Glueball Spectrum 

The gluon self-coupling that distinguishes QCD from the Abelian QED admits 
the existence of purely gluonic bound states, or glueballs. Indeed, glueballs are 
the only true states in quenched QCD! The good agreement of the quenched 
light-hadron spectrum with the experimental value was perhaps not surprising; 
Zweig's rule tells us that hadronic decays involving the annihilation of the 
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Figure 8: Data[E] for (a): vector mesons, (b): S-like baryons and (c): A-like baryons is 
plotted against the average value of the masses of the component quarks m q = (?B 9 ,i + 
m q 2)/2 (vector meson) and m q = i + rh 9 2 +?Ti 9 3)/3 (baryons). Squares and circles 



denote the N 
ofEqs. (Ell 



data at = 6.2 and f) = 6.0 respectively, and the lines correspond to the fits 
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Figure 9: The hadron spectrum using the Wilpqn and non-perturbatively improved fermion 
actions is shown against a 2 at fixed mps/trivtJ The scale is set from the string tension. 
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Figure 11: Construction of improved glueball operators, with the aim of increasing the 
overlap with the ground state. 



initial quarks are highly suppressed. Zweig's rule also suggests that glueball 
mixing with quark states should be similarly suppressed. Thus a quenched 
calculation of the glueball spectrum is very important, and can yield crucial 
information to guide experiment. 

The calculation of the glueball spectrum has been plagued by two inter- 
related problems. Firstly, the glueballs are relatively heavy and thus the cor- 
relation functions die rapidly at increasing temporal separations. Secondly the 
glueball correlators are subject to large fluctuations independent of separation. 
In consequence, the signal-to-noise ratio for glueball correlators is very poor. 

Calculations of the spectrum using the standard Wilson gluon action, 
Eq. (|TT|), have emphasised the construction of improved gluonic operators 
which more correctly describe the ground-state glueball wave function, as il- 
lustrated in Figure [Ll[ In the case of the determination of the spectrum of 
states at rest, the rotation group used in the construction of the continuum 
glueball operators is reduced to the cubic group of the lattice. The different 
components of, for example, the J PC = 2 ++ glueball lie in the E ++ and T 2 ++ 
representations of the cubic group. As the continuum limit is approached, the 
masses obtained from the different representations of the cubic group should 
become degenerate, signifying the restoration of rotation symmetry. Despite 
the discretisation errors of the standard Wilson gauge action being 0(a 2 ), the 
lightest glueball state is subject to much larger Q(a 2 ) discretisation errors than 
the spectrum for quark states, as shown in Figure |l2pil 

To improve upon these calculations, Morningstar and PeardonEH chose 
to employ an C(a 2 )-improved gluon action, having discretisation errors of 



24 




Figure 12: The continuum extrapolation of the masses of the J PC = ++ and 2++ states 
The different plotting symbols for the 2 ++ states correspond to the lattice operators T2 
(octagons) and E (diamonds). The lines represent linear continuum extrapolations in a 2 . 
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Figure-il3: The quenched glueball spectrum is shown after extrapolation to the continuum 
limitE3 
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0(a 4 );note that it is only necessary to consider operators of even dimension in 
the construction of the gluonic part of the action. Unfortunately, the relatively 
large mass of the glueball states in lattice units allows only a couple of time 
slices to be used in extracting the glueball masses. They then observed that 
one could employ a relatively coarse lattice in the spatial directions to produce 
a reasonable approximation to the glueball wave function, whilst employing 
a finer lattice in the temporal direction to enable the isolation of the ground 
state and excited state masses in each channelo 

The anisotropic lattice is implemented through the choice of different cou- 
pling constants for the space-space and space-time plaquettes in the gluonic 
action Eq. [ll]; this involves a non-trivial tuning—since the bare couplings are 
renormalised. The resulting glueball spectrumlp after extrapolation to the 
continuum limit, is shown in Figure 

2.4 Exotic Hadrons 

The search for hadrons with excited glue is one of the primary goals of the 
N* programme at CEBAF. There has been a flurry of recent activity looking 
for exotic states, and in particular exotic mesons, in the lattice community; 
lattice gauge theory aficionados generally try to comprehensively understand 
the mesonic sector before venturing into the realm of baryons. 

Exotics and Hybrids 

Within the quark model, the charge conjugation C and parity P of a meson 
are related to the spin S and orbital angular momentum L through 

P= (-1) L+1 ;C = (-l) i+s . 

States not conforming to these relations are called exotics, and examples are 
the states with 

J PC = l-+,0+-,2+-. 

An exotic can be formed in two ways. Firstly, as a quark-antiquark-glue state, 
which we call a hybrid. Secondly as a bound state of two quarks and two 
antiquarks. In this section I will discuss lattice studies of hybrid states. These 
studies have been given increased impetusJaY-experimental observations of 1 
resonance states in the region of 1.4GeVoc3 

Hybrid interpolating operators 

The usual interpolating operator for the pion is 

O v {x,t) = Tp(x,t) l5 iP(x,t). (44) 
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In the case of the hybrid state 1 + , a possible interpolating operator would be 
©!-+ & t) = i>{x, (x, t)1>{3, t) (45) 

where Fij is a lattice discretisation of the field-strength tensor constructed in 
Figure |§|. In practice to get any sort of signal for hybrid mesons, it is necessary 
to use interpolating operators in which the quark and antiquark are separated 
in space. The signals are inevitably much noisier than for the pseudoscalar and 
vector meson states. To see why, we note that the hybrid correlator Ch y b(i) 
decays exponentially with the mass of the hybrid, 

C hyh (t) ~ e ~ m '- bt . (46) 

In contrast, the correlator for the variance is that of the square of the interpo- 
lating operator 

C a2 (t) =X](|O hy b(x,t)| 2 |O hyb (0,0)| 2 ). (47) 

X 

Typically, |Ohyb| 2 is an interpolating operator for two pions, and therefore the 
signal-to-noise ratio with increasing temporal separations increases as 

S1 g nal ^ e -(m hyb -m,)t_ ^ 
noise 

Since the masses of the hybrids are relatively large, the signal quickly is lost 
in the noise. In the case of the glueball calculations, the situation is particu- 
larly severe since the square of the glueball operator has a non-zero vacuum 
expectation value, leading to the constant noise alluded to earlier. 

There have been recent calculations of the light-quark hybrid spectrum, 
and in particular, of the 1 h state, in both the quenched approximationp3ED 
and in full QCDeal23 These calculations all find a 1 h mass around 2 GeV, 
far larger than that of the experimental candidates. A possible resolution 
is that this resonance is actually a four-quark state; we will return to this 
interpretation in Section ^. 

2.5 The N* Spectrum 

The measurement of the excited nucleon spectrum reveals the full SU (3) na- 
ture of QCD, and is a critical part of the experimental programme at CEBAF. 
The observed N* spectrum is shown schematically in Figure ^ The nu- 
cleon 7V(938) has positive parity; its parity partner, with negative parity, is 
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Figure 14: Schematic showing the observed N* spectrum, labelled by J 



the 7V(1535). The usual nucleon interpolating operators employed in lattice 
calculations are 

N a = e ijk {u i C~ l sd?)u k a . (49) 

On forming the time-sliced correlator, both positive and negative parity states 
can contribute to the correlation function. However, we can perform a parity 
projection 

C{t) = ^(0 | N a {x, t)(l ± 7o) Q /3^(0) | 0) (50) 

X 

to project out the forward-propagating positive (+ sign) or negative (- sign) 
parity states, with states of the opposite parity propagating in the negative 
direction. On a periodic lattice, our correlator contains both the forward- 
and backward-propagating states, and we rely on a sufficiently long temporal 
extent to our lattice to delineate the two parities. 

Recently, there have been two lattice calculation of this mass splitting. The 
first Ej employs a highly improved fermion action, the -D x 34 action of Hamber 
and WuEJ Not only do the authors extract the mass of the J = l/2~j_state, but 
also find a signal for the J = 3/2 _ state. The second calculationcj employs 
domain- wall fermions; here the authors argue that, since the J = l/2 + and 
J = l/2~ state are degenerate in an unbroken chirally symmetric theory, the 
use of an action having a possessing exact chiral symmetry, even at non-zero 
lattice spacing, is crucial in correctly extracting the splitting. 
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Figure 15: The masses of the nucleon (N) and itp-negative-parity partner in lattice units 
(a -1 ~ 1.9 GeV using the p mass to set the scale)c3 
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Both calculations find mass splittings between the positive- and negative- 
parity states in accord with experiment, though with still substantial system- 
atic and statistical uncertainties. The masses of the J = 1/2+ (nucleon) and 
J = l/2~ states using DWF at a fixed value of the lattice spacing are shown 
in Figure [l5|. 

3 Hadron Structure 

As well as enabling the calculation of the hadron spectrum, lattice QCD en- 
ables the study of the distribution of the quarks and gluons within hadrons. 
Information about these is contained in the form factors, and in the quark and 
gluon structure functions. These calculations have in common the determina- 
tion of some (local) hadronic matrix element, so we will begin this section with 
a discussion of the lattice technology of determining matrix elements. 

3.1 Hadronic Matrix Elements 

The paradigm calculation is that of / w , the pion decay constant defined through 

(0\A fl \n)=ip ti U, (51) 

where = tpj^j^ij} is the axial vector current. This matrix element we can 
obtain as a by-product of the determination of m^, and our analysis will follow 
that of Section |2.2[ We construct the correlator 

C(i)=X>|4r (x,i)4 at (0)|0.), (52) 

X 

where A lat is a lattice discretisation of the continuum axial- vector current. 
Inserting a complete set of states between the two interpolating operators, we 
obtain 

= EE(^j3 / 2 ^e-^+^(0|4 at (0)|P(p))(P(p)|Ai att (0)|0) 

= E/ ^e^ Et S^Hp){0\Ar\P(p)){P(p)\AT f \0) 
*Z2? - m^ff 2 e~ m - * + excited states (53) 
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The lattice discretisation has provided both a cut-off, and a renormalisation 
scheme. We, of course, want / w is some familiar continuum scheme, such as 
MS. To provide us with that, we need to compute the matching coefficient Za 
for the axial vector current, such that 

A* = Z A A ut . (54) 

The determination of the matching coefficient Z is one of the most delicate, 
and generally onerous, tasks of any calculation of hadronic matrix elements. 
The lattice formulation reproduces continuum QCD as the lattice spacing ap- 
proaches zero, and thus the anomalous dimensions of the operators in the 
lattice formulation matches those of the continuum operators. However, the 
replacement of the continuum Lorentz symmetry by the hypercubic symmetry 
of the lattice, and the lack of chiral symmetry in most fermion formulations, 
complicates the calculation of the matching coefficients enormously, even for 
such a simple operator as the axial vector current. In particular, the lattice 
allows mixing with higher dimension operators, combined with appropriate 
powers of the lattice spacing a. An important element of the improvement 
programme is finding a combination of lattice operators such that matrix ele- 
ments are free of 0(a), or higher, discretisation errors. 

In principle, we can compute Z in perturbation theory. Perturbation the- 
ory using the bare lattice coupling g as an expansion parameter apparently 
fails, resulting in very large perturbative corrections. The bulk of these large 
corrections can be identified as lattice artifacts, the "tadpole" contributions 
arising from the 0(g 2 ) term in the expansion of the link variable, Eq. @. These 
terms can effectively be resiimmed through the expansion in terms of a renor- 
malised coupling constantly This prescription is precisely that used in the 
specification of the tadpole-improved clover coefficient of Eq. [}(]. 

An alternative route is to attempt to determine the matching coefficients, 
and improvement coefficients, non-perturbatively through the imposition of 
chiral Ward identities!!!} Indeed, this prescription enables, in principle, the 
elimination of all 0(a) discretisation effects from on-shell quantities, providing 
appropriate renormalisation conditions can be found. 

Both these routes require considerable effort, and uncertainty in the calcu- 
lation of matching and improvement coefficients is a major uncertainty in the 
calculation of hadronic matrix elements. Let me conclude this subsection by 



noting that the chiral fermion actions of Section 1.6 are automatically 0(a)- 
improved, and admit a smaller degree of operator mixing. This may be prove 
a substantial advantage for these formulations. 
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3.2 Nucleon Form Factors 



The electric and magnetic form factors of the nucleon are among the simplest 
quantities that contain information about the structure of the nucleon, and 
are measured in electron proton scattering. They are related to the matrix 
elements of the vector current J p through 



(p',s' | J M (g) \p,s) = u(p',s') 



^F 1 (q 2 )+ l C7^^—F 2 (q 2 ) 



u P {p,s), (55) 



where q = p — p' is the momentum of the photon probe. Note that -Fi and F% 
satisfy 

Fi(0) = l; f 2 (0) = A*-l 

where the former result expresses current conservation, and \x is the magnetic 
moment of the nucleon. For point-like particles, both quantities would be 
constant, and therefore they are a measure of the spatial extent of the nucleon. 
Rather than quoting these quantities directly, it is usual to form the Sach's 
form factors 

G E (q 2 ) = F 1 (q 2 ) + J^F^q 2 ) (56) 
G M (q 2 ) ^F 1 (q 2 )+F 2 (q 2 ), (57) 

where we note that q 2 is space-like. 

Phcnomcnologically, it is usual to parameterise the form factors through 
the vector dominance model by a dipole fit 

gW) - gim 2 )/^ ~ gUi 2 ) - - — X —-^2 

(1 - q 2 /m^) 

G n E (q 2 ) ~ 0. (58) 

A recent lattice determination of the proton electromagnetic form factors@ is 
shown in Figure 3.2, together with the experimental datau-3 The lines are fits 
to the data using the dipole forms of Eq. [58[ 

3.3 Hadronic Structure Functions 

The hadronic structure functions, describing the distribution of quarks and 
gluons inside, say, a nucleon in inclusive processes, are related to the hadronic 
tensor 

^ = h / ** e **( N fa a )\W J MW{p,a)), (59) 
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Figure 16: The solid points are lattice determinations of themfectric and magnetic form 
factors of the proton in tha— quenched approximation to QCDca The open circles are the 
experimental measurements Cj The lines represent the dipole fits of Eq. |58|to the experimental 
and lattice data. 
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where J ^ is the electroweak current, and p, s are the nucleon momentum and 
spin respectively. Decomposing W^ v according to the possible Lorentz struc- 
tures yields four structure functions, two spin-averaged, Fi2(x,Q 2 ) and two 
spin-dependent, gi.2(x,Q 2 ), where x is the Bjorken variable, and Q 2 = —q 2 . 
Phenomenologically, these are determined in Deep Inelastic Scattering, and 
parameterised at some reference energy scale. 

Near the light cone, x 2 ~ 0, the structure functions can be expanded using 
the operator product expansion (OPE) in terms of the matrix elements of 
certain local operators, of twist (dimension - spin) two, together with Wilson 
coefficients calculable in perturbation theory. Historically, it is these matrix 
elements that are computed in lattice simulations, since the non-zero lattice 



spacing in our calculations precludes the measurement of the currents in Eq. 59 
at sufficiently small separations, and more fundamentally it is unclear how to 
extract this quantity at light-like separations in Euclidean space. 

The local matrix elements are related to the ir-moments of the structure 
functions. In principle, we can recover the full x-dependence of the structure 
functions by measuring increasing moments. The simplest operators are the 
non-singlet operators for both the unpolarised and polarised structure func- 
tions, 

CV,...^ = irf^iD^ . ..iD^rip 

0jL. PB = ftftW D i* ■ ■ ■ iD^n/), (60) 

where symmetrisation of indices and removal of traces is understood, and the 
t is an SU(2) flavour matrix. The first faaLjapments moments for the nucleon 
have been measured by several groups e3'c3oH 

The lowest moment of the unpolarised quark distribution has a particularly 
simple interpretation in terms of the momentum fraction carried by the quarks. 
Figure ^ shows the first-moment of the u and d quark distributions in the 
quenched approximation^! both distributions are somewhat higher than phe- 
nomenological expectations. Unfortunately, it is unclear the extent to which 
the calculation of the first few moments of the structure functions enables a 
useful picture of the ^-dependence of the structure functions, and a means 
of performing a direct computation of the hadronic tensor, Eq. |5^, would be 
invaluable. 

The study of the flavour-singlet sector involves consideration of both the 
quark and gluon distributions, which mix. Computationally, these studies are 
much more demanding, but of tremendous interest since they venture beyond 
the simple valence picture of the nucleon. More recently, there have been 
attempts to investigate the role of higher-twist contributionspl to which up- 
coming experiments at JLAB should be particularly sensitive. 



35 




^ 1 1 1 1 1 1 1 1 1 1 1 

0.02 0.04 0.06 0.08 0.1 0.12 

am 

q 



Figure 17: The first moment of the unpolp-jsed u and d quark distributions in the nucleon 
as a function of the quark mass at /3 = 6.0Lj The circles and crosses denote results wjfcb the 
Wilson and clover fermion actions respectively. The bursts denote the CTEQ results 
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4 The Nucleon-Nucleon Interaction 



In the preceding, we have used lattice gauge theory to acquire a fundamen- 
tal understanding of the internal structure of an isolated hadron from first 
principles. It is natural to aim at a similar understanding of the interactions 
between hadrons, and in particular of the nucleon-nucleon interaction, the very 
foundation of nuclear physics. 

Understanding the strong interaction in multi-hadron systems from lat- 
tice QCD is a notoriously difficult problem. Multi-hadron states involve the 
computation of a four-point function and are relatively massive, and therefore 
the corresponding correlation functions quickly vanish into noise at increas- 
ing temporal separations. Furthermore, multi-hadron systems are large, and 
therefore the spatial extent of the lattice needs to be correspondingly larger 
than that used in hadronic spectroscopy. Finally, the use of a Euclidean lat- 
tice obscures the the extraction of the phase information of the full scattering 
matrixes Despite these difficulties, the problem is fundamental and compelling. 

Historically, there have been two approaches to this problem within lattice 
QCD. The first aimed to extract certain quantitative parameters, of the hadron- 
hadron interaction by direct lattice simulation. Luschc£3a£3 exploited the 
finite-size dependence to extract a discrete set of s-wave scatt^cjng lengths. The 
was thoroughly tested within an 0(4)-symmetric </> 4 jxic 
lengths have .been computed within QCD for pionsfiZrc 
Fiebig et al. explored the / = 2 tt — tt system by extracting a residual 
interaction potential, and then proceeded to compute the scattering phase 
shifts which were compared with experiment. 

The extraction of the s-wave scattering lengths of the tt — tt interaction has 
been encouraging, as illustrated in Figure |l8| where the s-wave scattering length 
for the isospin 1 = 2 pion-pion interaction is shownEj The investigation of the 
N — N system is more problematical. Here the scattering lengths are of the 
order of 10 fm, rather than less than 1 fm as is the case for the tt— tt interaction. 
Therefore, while lattice calculations do indeed find scattering lengths for the 
N — N interaction considerably larger than those for the tt — tt and tt — N 
interactional this approach is limited by our present inability to simulate on 
lattice sizes of the order of 10—20 fm, and at physical values of the pseudoscalar 
mass. 

The second approach is motivated by the realisation that important insight 
into the nucleon-nucleon interaction can be gleaned by studying the interac- 
tions of a much simpler system, that of two heavy-light mesons, with static 
heavy quarksEH Such a system exhibits most of the salient features of the 
nucleon-nucleon system, namely quark exchange, flavour exchange and colour 
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Figure 18: The 1 = 2 pion scattering length obtained in the quenched approximation to 
QCD is shown as a function of the lattice spacing, and in the continuum limitH Also shown 
as the open circle is the current-algebra prediction for m-^ = 140 MeV. 
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Figure 19: The connected (flavour-exchange) and disconnected contributions to the interac- 
tion between two heavy-light mesons are shown. The solid lines denote the heavy quarks, 
whilst the curly lines denote the propagation of the light-quarks. 
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polarization. The interaction between the mesons can be understood by the 
study of the four-point functions, shown in Figure |l9|: 

C^(x,r;y,s) = (0 1 (y)0 2 (s)d\(x)dl(r)). (61) 

Here each operator can either be that for a pseudoscalar (P) or for a vector 
(V). The range of the interaction can be assessed by forming the z- and i-sliced 
sum, 

C(R,T)= G^x^-y.s), (62) 

where 

x = (0,t 1 ),y=(d,t 1 +T) 

r = (aj_, R, t 2 ), s = (aj_, R, t 2 + T). 

At large distances R and times T, the correlation function is dominated by the 
lightest state | n), of mass M n , that can be exchanged between the mesons, 
with 

C(Je,T)~e-- Mi ' Jl <0|O 1 (5,r)6l(5,0)|n) x (n|O 2 (0,T)<5|(0,0)|0). (63) 

Can these correlation functions be constructed from the elemental quark prop- 
agators of Eq. p6| ? The connected diagram of Figure [l9] requires the evaluation 
of the propagator from one point on every time slice to every point on the 
lattice; this is manageable. Unfortunately, the computation of the discon- 
nected diagram requires the evaluation of all-to-all propagators. Nonetheless, 
the study of the connected diagram alone is valuable. It describes the flavour- 
exchange interaction, and lattice simulations have shown that the interaction 
is indeed mediated by meson exchange at large distances, and furthermore 
that the quantum numbers of the exchanged particle are in accord with naive 
expectations; for the process PP — > PP a vector meson is exchanged, whilst 
for the process PV — > VP a pseudoscalar particle is exchanged^ 

The especially attractive feature of heavy-light systems is that the heavy 
quarks admit the definition of a relative coordinate, and thereby a local adia^ 
batic potential. Recently, exploratory studies have been made of this potentially 
and evidence for nuclear binding sought, as illustrated in Figure The inves- 
tigation of this potential is important and more feasible than the measurements 
of the scattering lengths, because the large scattering lengths for the N — N 
system come from a short-range potential. Furthermore, by exploring the po- 
tential, we can discover the relative importance of gluon and meson exchange 
contributions at various distances. Understanding the nature of this poten- 
tial is also crucial to spectroscopy; are exotic mesons predominantly quark- 
antiquark-glue hybrid states, or four-quark states? 
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Figure 20: Adiabatic potential between two heavy-light mesons, with the light quarks in an 
isospin 1 = and spin S = configuration, from reference^. Measurements are obtained in 
the quenched approximation to QCD, on lattices at /3 = 5.7. The separations are measured 
in units of Ro = 0.53 fm, and the different plotting symbols correspond to different lattice 
sizes. 
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5 Conclusions 



In these lectures I have tried to convey the power of lattice gauge calculations 
as a means of understanding hadronic physics. There are very many areas that 
I have not addressed - topology and the role of instantons, finite-temperature 
and finite-density phase transitions, and, reaching beyond hadronic physics, the 
calculation of weak-interaction matrix elements and even supersymmetry and 
gravity. I trust that I have convinced you that lattice gauge theory provides 
not only an ab initio tool for obtaining quantitative results (the spectrum, 
form factors etc.), but also provides a means of increasing our conceptual 
understanding of the strong interaction, for example through the study of the 
nuclear- nuclear force. 
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